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Summary  of  Results  and  Accomplishments 


Technical  summary. 

Electro-optic  polymer  relaxation 
Objective: 

To  better  understand  the  orientational  decay  mechanisms  of  poled  second-order  NLO 
polymers  and  to  develop  a  predictive  description  of  the  long  time  scale  (in  terms  of  tens 
of  years)  orientational  decay  of  poled  polymer  systems. 

Approach: 

An  electric  field  poling  process  is  used  to  create  noncentrosymmetric  order  in  polymer 
hosts  containing  dipolar  nonlinear  optical  chromophores.  The  resulting  materials  have 
been  found  useful  for  producing  a  variety  of  nonlinear  optical  and  electro-optic  devices. 
The  easy  processability  into  thin  films  make  such  material  ideal  for  waveguide  electro¬ 
optic  devices.  A  number  of  companies  are  now  developing  such  commercial  devices. 

A  major  issue  in  these  materials  arises  from  the  nonequilibrium  nature  of  the  polymers 
resulting  in  the  eventual  decay  of  the  nonlinear  optical  response,  which  is  due  to  the 
reorientation  of  the  aligned  chromophores  into  a  randomly  oriented  state.  Therefore, 
much  research  has  centered  on  understanding  the  decay  mechanism  with  the  goal  to 
guide  the  design  of  new  and  more  stable  materials.  To  this  end,  our  program  aims  to 
develop  a  predictive  model,  which  is  tied  to  relevant  material  parameters  such  as  the 
glass  transition  temperature,  and  the  relaxation  behavior  at  this  temperature.  As  materials 
improve,  and  stability  increases,  the  need  for  accelerated  testing  increases.  Thus,  we  seek 
to  develop  models  and  measurement  regimes  that  probe  behavior  at  high  temperatures 
(thus  short  times)  and  can  be  connected  theoretically  to  low  temperature  (thus  long  time) 
behavior. 

The  orientational  relaxation  of  nonlinear  optical  chromophores,  which  are  aligned  within 
a  polymer  using  electric  field  poling  techniques,  is  highly  coupled  to  the  relaxation 
behavior  of  the  host  polymer  glass.  Polymer  glasses  are  well  known  to  display  highly 
dispersive  relaxation  behavior  as  the  glass  tends  toward  equilibrium.  This  dispersive 
behavior  is  manifested  in  relaxation  processes  being  spread  over  all  time  scales.  Indeed, 
it  is  believed  that  such  behavior  is  self-similar  over  a  very  wide-range  of  time  scales. 
Additionally,  it  has  been  found  that  time  and  temperature  behavior  can  be  superimposed. 
We  wanted  to  develop  measurement  techniques  that  probe  the  relaxing  nonlinear  optical 
chromophore  directly,  and  to  be  able  to  probe  the  coupling  to  the  polymer  over  as  wide  a 
time  and  temperature  scale  as  possible.  By  doing  this,  we  could  relate  the  chromophore 
relaxation  to  the  dielectric  relaxation  in  the  polymer.  The  concept  of  superposition  of 
time  and  temperature  will  allow  for  the  development  of  accelerated  lifetime  testing  when 
an  appropriate  model  is  developed. 
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We  studied  the  relaxation  behavior  over  many  decades  of  time  using  frequency  and  time 
domain  techniques.  Based  on  these  results,  we  were  better  able  to  probe  a  wide 
temperature  range  to  ascertain  the  characteristic  times  and  dispersion  parameters  of  the 
relaxation  phenomena.  We  were  then  able  to  determine  which  models  best  fit  the  data, 
and  then  to  use  these  models  to  predict  behavior  at  different  temperatures  and  time  scales. 


Results: 

We  have  assembled  the  temperature  control  system  shown  in  Figure  1  which  can  access 
temperature  from  liquid  nitrogen  temperature  to  about  1 50  C.  We  have  a  separate  cell 
which  can  access  temperatures  to  300  C.  The  apparatus  in  Figure  1  consists  of  a  liquid 
nitrogen  dewar  with  a  resistance  heater  which  is  temperature  controlled.  This  system  has 
been  built  to  be  capable  of  carrying  out  electrical  measurements  concurrent  with  the 
nonlinear  optical  measurements.  Thus,  we  have  the  capability  to  measure  the  dielectric 
constant,  resistivity,  thermally  stimulated  depolarization,  and  isothermal  depolarization. 
These  measurements  will  be  able  to  characterize  the  relaxation  behavior  of  the  polymer 
over  an  extremely  wide  range  of  times  so  as  to  best  determine  the  appropriate  model  (see 
above)  which  describes  the  relaxation. 


Figure  1.  Dewar  system. 
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We  built  a  chielectric  spectroscopy  apparatus  for  measuring  the  orientation  response  of 
the  nonlinear  optical  chromophore  directly  in  the  frequency  domain.  This  is  a  new 
technique,  which  we  have  improved  by  using  a  mode-locked  CW  laser  and  a  lock-in 
amplifier  to  measure  the  frequency  response  of  electric  field  induced  second  harmonic 
generation  over  a  broad  frequency  range.  The  mode-locked  Ti: Sapphire  laser  is 
described  Richard  Dureiko’s  Ph.D.  Thesis.  The  experimental  setup  is  shown  if  Figure  2. 
The  entire  chielectric  experiment  is  described  in  detail  in  Appendix  1 .  This  technique 
was  found  useful  for  measurements  near  and  above  the  glass  transition  temperature. 


lock-in 

int 

osc. 

amplifier 

in 

10  mHZ-  40  KHz 


power  amplifier 


Figure  2.  Chielectric  spectroscopy  setup. 


We  measured  the  dielectric  and  chielectric  response  near  and  above  the  glass  transition  in 
the  frequency  domain  for  a  series  of  doped  and  covalent  polymers  based  on  methyl 
methacrylate  and  in  cross-linked  materials  provided  by  L.  Dalton  (see  R.  Dureiko  thesis). 
We  also  measured  the  dielectric  and  chielectric  response  in  the  time  domain  below  the 
glass  transition  temperature.  Using  all  four  techniques,  we  were  able  to  probe  the 
response  on  many  time  scales  and  temperature  regimes. 

We  analyzed  the  isothermal  data  in  terms  of  a  three  parameter  relaxation  model  as 
adapted  from  the  model  of  Dissado  and  Hill.  This  analysis  provides  the  dispersion 
parameters  and  characteristic  times  of  the  relaxation.  The  data  was  collected  at  a  range  of 
temperatures  both  above  and  below  the  glass  transition  temperature,  and  from  it  we  could 
conclude  that  a  three-parameter  model  including  dispersion  and  asymmetry  was 
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Temperature  Dependence  of  Relaxation  Times 


Figure  3.  Reduced  temperature  dependence  plot  for  x  scaled  to  xg  for  the 
three  polymer  systems.  The  solid  line  represents  the  Adam-Gibbs  model 
based  on  enthalpy  measurements  from  DR1-MMA.  The  dashed  line 
represents  an  Arrhenius  fit  based  on  data  points  below  the  glass 
transition,  while  the  dotted  line  represents  Fulcher-Tammann-Hesse  or 
equivalently  WLF  behavior  above  the  glass  transition. 

necessary  to  fit  the  data.  A  single  constant  diffusive  model  and  even  the  widely  used 
KWW  function  was  inadequate  to  fit  the  data  over  the  many  time  decades  we  were  able 
to  measure.  Thus  our  program  of  time  and  frequency  domain  measurements  scanning 
many  decades  of  time  was  able  to  differentiate  between  models  for  the  isothermal 
relaxation. 

Finally,  we  applied  the  Adam-Gibbs  model  to  the  isothermal  data  at  many  temperatures. 
Again,  our  technique  provided  for  measurement  over  many  time  scales  which  also  allows 
determination  of  the  characteristic  times  over  a  wide  temperature  range.  This  allowed  for 
good  comparison  to  the  Adam-Gibbs  model.  In  this  model,  the  temperature  behavior  is 
well-described  both  above  and  below  the  glass  transition.  The  model  agrees  with  other 
models  which  are  only  applicable  over  limited  temperature  ranges.  Thus,  the  model  is 
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able  to  connect  high  temperature  (short  time)  measurements  with  low  temperature  (long 
time)  measurements  and  provides  the  basis  for  a  predictive  model.  In  addition,  the 
characteristic  times  were  normalized  to  that  at  the  glass  transition  and  the  temperature  to 
the  glass  transition  in  a  data  plot  as  shown  in  Figure  3.  Using  this  technique,  predictive 
capability  results,  and  further,  it  is  seen  that  similar  polymers  fall  onto  a  single  curve  for 
quasi-universal  behavior.  Thus  predictions  can  be  made  on  closely  related  systems  using 
the  glass  transition  as  a  scaling  parameter.  This  allows  for  guidance  in  new  material 
development.  Details  in  Appendix  1. 

We  believe  that  we  have  attained  the  major  goals  of  the  program.  We  have  identified  a 
model  and  measurement  regime  appropriate  for  accelerated  failure  tests,  which  can  aid  in 
the  development  of  more  stable  electro-optic  polymers.  In  doing  this,  we  have  developed 
a  new  measurement  technique,  chielectric  spectroscopy,  which  can  be  applied  to  study 
basic  questions  in  polymer  physics.  Since  the  technique  is  sensitive  only  to  the 
chromophores,  it  can  be  used  to  study  whether  the  relaxation  is  a  collective  effect  of  the 
dipoles. 

Laser  system  development 

We  built  a  new  tunable  laser  source  for  harmonic  generation  experiment.  The  system 
consists  of  an  active-passive  modelocked  Nd:YAG  laser  with  negative  feedback  in  the 
oscillator  to  produce  a  train  of  about  50  20  picosecond  pulses  during  each  lampflash.  The 
laser  contains  a  double  pass  amplifier  followed  by  frequency  tripling  crystals  whose 
output  pumps  a  parametric  oscillator.  It  is  the  first  subpicosecond  Nd:YAG  based  laser 
source.  It  uses  a  synchronously  pumped  optical  parametric  oscillator  (SPOPO)  which 
when  properly  tuned  generates  subpicosecond  tunable  laser  light.  Pulses  as  short  as  400 
fs  have  been  observed.  The  system  uses  BBO  crystals  and  is  tunable  from  450-2500  nm. 


Figure  4.  Autocorrelation  trace  of  subpicosecond  laser  pulses  from  the 
SPOPO  depicting  pulse  compression. 
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Short  pulses  are  generated  using  group  velocity  walk  off  in  the  slightly  detuned  SPOPO. 
We  have  modeled  the  pulse  generation  mechanism  using  the  coupled  wave  equations  in 
the  nonlinear  optical  crystal  and  have  studied  the  mechanism  responsible.  The  pulse 
compression  involves  the  group  velocity  walk  off  and  the  depleted  parametric 
conversion.  It  is  believed  that  the  compressed  pulses  are  3  wave  temporal  solitons. 

This  unique  source  was  utilized  to  measure  the  dispersion  in  third  harmonic  generation  in 
squaraine  dyes.  In  addition,  we  have  recently  used  this  laser  as  the  core  of  the  only 
widely  tunable  hyper-Rayleigh  scattering  apparatus.  We  have  used  it  to  study  the 
dispersion  of  hyper-Rayleigh  scattering,  and  have  identified  dephasing  as  a  problem  in 
the  measurement  of  near-resonant  hyper-Rayleigh  scattering.  Most  recently,  we  have 
measured  Kleinman-symmetry-breaking  in  organic  molecules,  which  are  candidates  for  a 
new  class  of  “quadrupoled”  polymeric  materials. 

Crosslinked  polyimide  materials  for  integrated  optics 

We  have  been  studying  a  promising  new  class  of  materials  based  on  crosslinkable 
polymers  developed  by  Amoco  Chemical  Co.  The  crosslinking  makes  the  materials 
desirable  for  stable  electro-optic  films,  and  their  photosensitivity  makes  device 
processing  extremely  attractive. 

We  studied  the  structure  property  relationships  associated  with  waveguide  losses  and 
successfully  developed  very  low  loss  materials  which  are  now  being  investigated  by  a 
number  of  companies  for  waveguide  passive  interconnect  applications.  We  have 
successfully  poled  nonlinear  optical  molecules  in  these  polymers  and  have  observed 
stable  electro-optic  properties. 

We  have  investigated  poling  through  multilayer  waveguide  structures  and  have 
developed  a  technique  to  dope  the  cladding,  which  concentrates  the  poling  field  in  the 
active  layer.  This  approach  led  to  nearly  a  tenfold  increase  in  electro-optic  coefficient. 

Devices  in  similar  materials  are  now  under  development  by  Lockheed-Martin.  Details  are 
given  in  Appendix  2. 
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Abstract 


Chielectric  spectroscopy,  which  is  frequency  domain  electric  field  induced  second  harmonic 
generation,  was  developed  to  study  the  nonlinear  optical  response  near  the  glass  transition  temperature  at 
short  time  scales.  These  measurement,,  along  with  time  domain  second  harmonic  generatton  and  time  and 
frequency  domain  dielectric  measurements  were  used  to  characterize  the  decay  of  poling-mduced  electro- 
optic  properties  of  guest-host  and  side-chain  methacrylate  polymers  having  glass  rransirion  temperatures  in 
the  range  of  90<  T,  <  1 25  'C.  Time  and  frequency  dependent  data  were  modeled  with  the  model  of  Dtssado 
and  Hill.  The  temperature  dependence  of  the  time  constants  from  rhe  polymeric  systems  were  then 
compared  and  a  scaling  model  for  predicting  useful  lifetimes  of  poled  eleetro-optic  media  is  discussed. 
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1.  Introduction 


Over  the  pas.  decade,  a  great  deal  of  research  has  been  focused  on  the  development  of  nonlinear 
optical  (NLOI  and  electro-optic  (EO)  polymeric  materials  because  of  their  potential  in  applications  ranging 
from  second  harmonic  generation  (SHG)  and  phase  conjugation,  to  optical  switching  and  electro-optic 
modulation.1,2'3  Due  to  their  large  nonlinear  optical  properties,  ease  of  processing,  and  expected  low  cost, 

polymers  are  leading  candidates  for  electro-optic  materials. 

Electric-field  poling  has  become  a  favored  technique  for  producing  noncentrosymmetric  polymeric 
electro-optic  materials.  Here,  nonlinear  optical  molecules  incorporated  in  a  polymer  host  are  aligned  in  an 
electric  field  when  in  a  highly  mobile  state,  and  the  resulting  polar  orientational  order  is  “locked-in”  by 
cooling  or  chemical  reaction.  The  decay  of  the  electro-optic  and  nonlinear  optical  properiies  arises  from 
the  relaxation  of  the  polar  orientational  alignment  originally  induced  in  the  electric  field  poling  process. 
The  polar  orientaiionai  order  eventually  reverts  to  the  original  isotropic  alignment  present  prior  to  poling, 
and  is  due  to  the  orientational  mobility  of  the  aligned  nonlinear  optical  chromophores  as  well  as  to  the 
relaxational  motion  of  the  polymer  host  in  its  non-equilibrium  state.  In  order  for  a  polymeric  device  to  be 
viable,  i,  must  retain  its  NLO  properiies  over  the  course  of  the  required  lifetime  of  the  device,  which  may 
be  several  years.  The  usefulness  of  devices  based  on  these  poled  polymeric  materials  is  dependent  upon  the 
temporal  and  temperature  stability  of  the  NLO  properties  of  the  material. 

Since  the  earliest  poling  experiments,  the  instability  of  the  orientational  order  was  recognized,  as  were 
paths  to  the  development  of  new  materials  which  would  enhance  the  stability.  A  great  deal  of  effort  has 
gone  into  developing  new  stable  materials  as  well  as  into  studies  aimed  at  understanding  the  process 
underlying  the  relaxation.  It  has  been  found  that  a  close  relationship  exists  between  the  decay  of  the 
electro-optic  coefficient  and  other  relaxation  processes  in  polymers,  and  so  has  been  able  to  draw  on  that 
literature.  A  number  of  approaches  have  been  taken  to  describe  the  process,  and  great  progress  has  been 
made  in  developing  a  predictive  model  for  the  lifetime  of  electro-optic  polymers.1 

In  this  paper,  our  aim  is  to  move  closer  to  a  universal  predictive  model  and  a  regime  of  measurements 
useful  in  predicting  lifetime  from  measurements  taken  in  a  short  time.  We  have  approached  the  subject 
with  the  goal  of  unifying  the  disparate  approaches  that  have  been  used  in  the  past  by  seeking  models  for  the 
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process  which  subsume  previously  published  models.  We  begin  with  a  review  of  the  literature  in  the 
context  of  models  and  experimental  techniques  and  then  describe  these  in  relation  to  the  somewhat  more 
general  models  we  propose  here.  In  addition,  we  describe  a  recently  developed  experimental  technique, 
“chielectric  spectroscopy",  which  measures  electric  field  induced  second  harmonic  generation  in  the 
frequency  (modulating)  domain.  This  technique  is  analogous  to  dielectric  spectroscopy  except  that  it 
directly  probes  only  the  nonlinear  optical  chromophore.  With  this  technique  and  time-domain  second 
harmonic  generation  decay  measurements,  a  wide  frequency  range  of  the  orientational  response  can  be 
probed,  corresponding  to  a  wide  range  of  temperatures  and  times.  Finally,  we  describe  measurements  on 
acrylate  polymers  in  the  context  of  our  proposed  experimental  regime  and  models  with  the  goal  of  a 

predictive  model. 


2.  Orientation  Relaxation  Models 

As  mentioned  previously,  a  requirement  for  second-order  nonlinear  and  Pockets  materials  is  that  the 
material  lacks  an  inversion  center  (noncentrosymmetric).  The  electric  field  poling  process  has  been 
developed  to  impose  this  symmetry  on  otherwise  isotropic  polymer  materials  containing  nonlinear  optical 
chromophores.  In  this  process,  polar  order  is  imposed  which  can  be  described  by  an  oriented  gas  model, 
for  which  the  Pockels  coefficient  is  given  by,5 

rx{-co:co, 0)  =  A^cos3  ^/"/"/"^(-tutro.O)  (0 

where  N  is  the  number  density  of  nonlinear  optical  chromophores.  the  fs  are  local  field  factors  at  the 
specified  frequencies,  /?.=  (-<y;ru,0)  is  the  molecular  hyperpolarizability,  and  <cos30>  is  a  thermal  average 
which  depends  on  the  angle  between  the  molecular  dipole  and  the  poling  field  direction,  the  strength  of  the 
molecular  dipole  moment,  p,  and  the  strength  of  the  applied  field.  E.  This  thermal  average  is  the  order 
parameter,  and  decays  with  time  due  to  the  orientational  mobility  of  the  chromophore  and  the  relaxationai 
motion  of  the  polymer  host.  Once  the  field  is  removed,  the  chromophores  find  themselves  in  a  thermal 
non-equilibrium  environment,  and  will  tend  to  randomly  orient  themselves  with  respect  to  each  other. 
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resulting  in  a  decay  of  r„  over  time.6'7  This  relaxation  of  the  chromophore  dipoles  is  hindered  by  the 
polymer  matrix  which  surrounds  them. 

When  the  poling  field  is  turned  off.  one  can  detect  the  orientational  relaxation  through  a 
measurement  of  the  second  harmonic  or  electro-optic  coefficient.  A  description  of  the  decay  of  the  optical 
nonlinearity  has  two  parts  which  ultimately  provide  for  the  time  and  temperature  dependence  of  the  decay 
of  the  nonlinear  optical  properties.  In  the  first  part,  the  time  dependence  of  the  isothermal  decay  is 
modeled  using  chromophore  reorientation  models,  a  number  of  which  are  described  below.  These  models 
require  experimentally  determined  parameters  which  characterize  the  decay  function.  In  polymer  systems, 
two  or  three  parameters  should  suffice,  although  some  models  contain  more.  Generally,  one  of  the 
parameters  describes  the  characteristic  time  of  the  system,  while  the  other  one  or  two  parameters  describe 
the  dispersive  nature  of  the  decay  exhibited  in  viscoelastic  systems  like  polymers.  These  parameters  are 
characteristic  of  the  materials  and  the  processing  history.  It  is  well  known  from  the  dielectric  literature, 
that  the  time  dependence  can  be  probed  in  the  frequency  domain.  We  further  develop  this  approach  here 
with  the  measurement  of  chielectric  spectroscopy  where  frequency  dependent  electric  field  induced  second 
harmonic  generation  is  studied.  Measurements  in  the  time  and  frequency  domains  are  equivalent  and  are 
related  through  appropriate  mathematical  transforms.  In  the  second  part  of  developing  a  model,  a  physical 
picture  of  the  material  properties  contributing  to  the  decay  is  drawn.  This  part,  the  relaxation  model,  will 
describe  the  temperature  dependence  of  the  parameters  describing  the  chromophore  reorientation  which  are 
obtained  through  a  fit  to  the  isothermal  time-dependent  data.  Here  the  statistical  physics  is  introduced  and 
related  to  the  temperature  dependence  of  the  chromophore  reorientation  parameters.  We  will  now  review 
the  models  used  in  the  literature  and  introduce  our  models  of  choice  here,  namely  the  Dissado-Hill  model 
for  the  chromophore  reorientation,  and  the  Adams-Gibbs  model  of  the  relaxation. 

2.1.  Chromophore  reorientation 

2.1.1.  Previously  studied  models 

One  of  the  simplest  decay  models  is  the  Debye  model.8  In  this  model,  chromophore  orientational 
relaxation  is  described  by  a  single  exponential  function. 
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(2) 


P(t)=  Pn  exp(-//r) 

where  P  represents  the  nonlinear  optical  polarization  and  t  is  a  characteristic  relaxation  time.  Thus,  the 
square  root  of  the  second  harmonic  intensity  versus  time  could  be  fit  to  this  function  where  applicable.  This 
one  parameter  model  is  helpful  in  describing  relaxation  in  gases  and  dilute  low  viscosity  liquids.  The 
shortfall  of  this  simple  model  is  borne  out  in  experimental  data  that  clearly  indicates  that  it  fails  to  describe 
the  time  dependence  of  the  orientational  decay  in  viscoelastic  materials,  as  has  been  pointed  out  by  a 
number  of  workers.',l<uu2  The  Debye  model  fails  to  capture  the  dispersive  nature  of  the  relaxation  in  such 
materials.  The  failure  of  this  model  can  be  understood,  for  example,  by  considering  that  the  individual 
chromophores  undergoing  reorientation  are  surrounded  by  a  complex  polymer  network  that  hinders  their 
rotations.  Due  to  this  network,  individual  chromophores  will  experience  different  molecular  environments 
that  hinder  their  relaxation,  which  means  that  one  should  expect  a  range  of  characteristic  relaxation  time 
constants  instead  of  just  one. 

A  second  model  addresses  this  aspect  by  taking  into  account  that  more  than  one  characteristic  time 
constant  is  needed  to  describe  the  orientational  relaxation  of  the  chromophore  molecules  in  the  polymer 
matrix.  Thus,  the  distribution  of  relaxation  times  (DRT)  improves  upon  the  simple  Debye  model.13' 

17- 18  This  model  is  similar  to  the  Debye  model  in  that  it  assumes  the  relaxing  dye  molecules  are  completely 
independent  of  each  other,  and  thus,  there  are  no  dipole-dipole  interactions.  It  differs  in  that  the  nonlinear 
optical  polarization  is  now  written  as  a  superposition  of  Debye  relaxations. 

P{t)=fjP,exp(-t/Ti)  (3) 

/=0 

where  the  t,'s  represent  the  individual  characteristic  relaxation  times.  This  can  also  be  written  for  a 
continuous  distribution  described  by  a  function  p{ r)  using, 

/>(/)  =  jp(r)exp(-//r)r/r  .  (4) 

0 

Another  model,  which  can  be  viewed  as  a  limiting  form  of  the  discrete  DRT  model,  is  the  bi- 
exponential  model  which  describes  the  decay  process  in  terms  of  only  two  characteristic  time  constants,  t,, 

and  t2,I92° 
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P(t )  =  P,  exp(- 1/ r, )  +  P;  exp(-  t/r2 )  (5) 

Researchers,  such  as  Goodson  er  a/.2'  and  Hampsch  er  al.”  have  found  this  type  of  model  fits  their  second 
harmonic  (SH)  decay  data  well  in  PMMA  polymer  systems  doped  with  DANS. 

Currently,  the  most  widely  used  empirical  decay  model  in  the  literature  is  the  Kohlrausch-Wilhams- 

Watts  (KWW)  or  ‘stretched  exponential’  model. 


(6) 


P(t)  =  P0  exp[-  {t/r Y  ] 

where  x  is  the  characteristic  relaxation  time  constant  and  P  is  a  parameter  that  varies  between  zero  and 
unity.  This  model  was  introduced  by  Kohlrausch  in  1847  to  describe  mechanical  creep  in  glassy  fibers23 
and  then  again  over  100  years  later  in  1970  by  Williams  and  Watts24  to  describe  dielectric  relaxation.  Since 
then,  this  model  has  been  used  by  several  researchers  to  interpret  SH  decay  data.  2”s'26'27  Although 
originally  applied  as  an  empirical  law.  researchers  such  as  Klafter  and  Shlesinger,28  and  even  more  recently 
Weron29  have  tried  to  bring  some  physical  meaning  to  this  law.  The  KWW  model  can  be  related  to  the 

DRT  model  in  that  we  can  write  it  as. 


exp(-//r)'a  =  jexp(-f/r)  p{r)  dr  0) 

o 

where  p(t)  is  .he  dis.ribucion  of  Debye  relaxation  times  introduced  above.”  In  this  interpretation,  p  is  a 
shape  parameter  of  the  distribution  function.  When  P  is  unity,  .he  KWW  function  reduces  to  the  Debye 
model.  Several  researchers  have  found  that  (1  increases  with  increasing  temperature,  but  no  good  model  of 
its  temperature  dependence  exists. “6‘ 31 

The  time  scales  over  which  experiments  have  monitored  chromophore  reorientation  range  from 
microseconds  in  DR1  doped  PS32  all  the  way  out  to  a  year  in  a  polyimide  side-chain  polymer,33  and  m  both 

studies,  the  KWW  was  used  to  tit  the  SH  decay  data. 

Another  empirical  model  put  forth  by  Van  der  Vorst  et  al. 20  is  a  modified  version  of  the  classical 

Debye  model.  In  this  model  the  characteristic  relaxation  time  constant  now  becomes  a  function  of  time, 

P(t)  =  P0  exp(-  r{t)  =r0+AtB  (8) 
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where  A.  B  and  x0  are  fit  parameters.  Adding  these  additional  parameters  to  the  model  allowed  a  better  fit 
to  their  experimental  data  on  DANS  side-chained  polymers,  however,  the  physical  justification  for  doing  so 
is  unclear. 

Several  researchers  have  put  forth  relaxation  models  based  on  chromophore  rotational  Brownian 
motion  in  the  polymer  matrix.33' 35  Liu  et  al.36  have  developed  a  Brownian  motion  model  that  also  accounts 
for  electric  field  effects,  and  polymer  restrictions  to  describe  their  corona  poled  DANS  doped  PMMA 
systems.  In  their  model, 

X~  ~  A'  f,  M  f:  (2  u)p:::  {cos3  0  p(9,f,  a ,  b , 90 ,  D)  sin  6  d6  (9) 

0 

where  f  is  the  local  field  factor.  N  is  the  number  of  rotating  chromophores.  and  is  the  second-  order 
polarizability  of  each  chromophore.  The  effects  due  to  the  electric  field  are  accounted  for  in  the  parameters 
a  and  b,  and  polymer  restrictions  are  represented  by  the  wedge  angle  90.  The  only  adjustable  parameter  in 
this  model  is  the  diffusion  constant,  D,  the  other  parameters  are  determined  in  independent  experiments.  It 
has  not  been  established  that  this  model  describes  the  dispersive  nature  of  the  relaxation. 

As  another  example  of  the  use  of  distributions  of  relaxation  times  (DRT)  ,  Verbeist  et  al.  have 
shown  that  the  Wagner  or  lognormal  distribution  function  is  useful  in  describing  chromophore  relaxation. 
In  this  model  the  decay  function  is  described  by  a  continuous  distribution  of  exponential  relaxation  times, 

Pit )  =  —  P  exp(-  u 2 )  exp[-  x exp(-  f3a  u)\du  (10) 

where  p0  is  the  half-width  of  the  Gaussian  distribution  at  the  1/e  point,  x  is  defined  as  x/xm,  with  xm  being 
the  relaxation  time  associated  with  the  maximum  of  the  distribution  function,  and  u  defined  as  ln(T/xm). 


2.1.2.  Dissado-Hill  Many  Body  Model 

Departing  from  the  empirical  models  already  reviewed  above,  we  now  briefly  outline  the  Dissado-Hill 
(D-H)  Many  Body  Universal  Relaxation  Model.j8  A  detailed  account  of  this  theory  and  the  experimental 

39 

evidence  that  supports  it  as  being  a  universal  model  for  dipolar  relaxation  has  been  made  by  Jonscher. 
This  model  is  on  sounder  footing  than  the  other  dispersive  models  above  in  that  it  is  capable  of  predicting 
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the  temperature  and  time  dependence  of  the  amplitude  of  the  response.  Although  this  model  was 
introduced  to  describe  dielectric  loss,  we  have  used  it  to  describe  chromophore  orientational  relaxation  in 
both  the  time  and  frequency  domain,  as  will  be  discussed  later  in  this  paper.  Unlike  other  decay  models 
which  assume  the  relaxing  dipoles  are  non-interacting,  the  D-H  model  assumes  just  the  opposite,  that 
dipole  relaxation  is  a  highly  cooperative  process.39  Additionally,  the  model  statistically  combines  behavior 
in  three  frequency  regimes.  In  high  and  low  frequency  regime,  the  dispersive  behavior  is  described  by  a 
simple  power  law.  Such  a  law  has  a  firm  basis  in  that,  unlike  the  other  empirical  models,  is  consistent  with 
Kramers-Kronig  analysis,  and  thus,  causality.  In  the  transitional  regime,  the  Dissado-Hill  model  yields  the 
temperature  and  frequency  dependence  of  the  response  amplitude  which  are  absent  in  other  models.  The 
assumptions  made  in  this  model  are  that  there  are  two  preferred  orientational  states  the  dipoles,  thus  giving 
rise  to  a  two  level  model.  The  potential  energy  diagram  is  shown  in  Figure  1. 

In  Figure  1.  A  is  the  height  of  the  energy  barrier  that  separates  the  two  preferred  states,  and  £,  denotes 
the  width  of  correlated  states.  When  an  external  electric  field  is  applied  to  the  dipolar  dielectric  media,  the 
two  potential  well  bottoms  are  split  by  an  amount  that  depends  not  only  on  the  strength  of  the  applied 
electric  field,  but  also  on  the  occupancies  of  the  two  preferred  energy  states.  There  are  three  types  of 
transitions  that  can  occur,  large  thermal  transitions,  small  flip  transitions,  and  finally  flip-flop  or  fluctuation 
transitions.  The  large  transitions  are  thermally  activated  Debye-like  transitions  and  thermally  assisted 
tunneling  processes  that  take  a  particle  in  one  of  the  preferred  states  across  or  over  the  potential  barrier  into 
the  other  preferred  state.  Small  flip  transitions  are  tunneling  transitions  between  the  two  preferred 
configurations  that  result  in  a  net  change  of  the  total  dipole  moment,  similar  to  the  type  of  transitions  Ngai 
e,  ai"  had  proposed.  The  final  type  of  transitions,  the  flip-flop  transitions  can  be  described  as  local  dipole 
moment  fluctuations  which  involve  the  simultaneous  tunneling  of  two  particles  between  the  two  preferred 
states,  but  in  opposite  directions,  resulting  in  no  net  change  in  the  total  dipole  moment. 

These  three  types  of  transitions  are  in  competition  with  each  other,  and  each  has  a  time  scale  in  which 
it  is  the  dominant  process.  At  times  much  shorter  than  l/cop,  where  cop  is  a  characteristic  frequency  in  the 
Dissado-Hill  model,  the  flip  transitions  are  the  dominant  processes.  At  times  approximately  1/  <o„  the 
dominant  processes  are  the  large  thermal  Debye-like  transitions.  Finally,  at  times  much  longer  than  1/  cop, 
the  flip-flop  or  fluctuation  transitions  take  over.  The  degree  to  which  the  flip  transitions  are  correlated  in  a 
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particular  dipolar  media  determines  the  shape  of  the  frequency  dependent  dielectric  loss  curve  at 
frequencies  higher  than  cy  Similarly,  the  degree  of  correlation  ot  the  flip-flop  transitions  determines  the 
shape  of  the  loss  curve  at  frequencies  below  o)p. 

The  complete  description  of  the  dielectric  loss  curve  as  a  function  of  frequency  is  a  combination  of 
these  three  processes,  and  can  be  written  in  terms  of  the  Gaussian  hypergeometric  function.  2Fi, 

z{a)=z(0t'  +  i6>/«>p),~l2Fk-n  ’  2~m  :  2-«;(l  +  ia>/a>pY) 


X"(co)oza)"-\  co»cop  (11) 

%"{co)  oc  Im(l  +  ico/cop ) 1 ,  co  =  cop 
Z"{co)occo"'.  C0«0)p 

Here  n  is  a  parameter  that  describes  the  degree  of  correlation  of  flip  transitions,  and  m  is  a  parameter 
that  describes  the  degree  of  correlation  of  flip-flop  transitions.  Both  parameters  can  range  from  zero  to 
unity.  As  noted  by  Jonscher59  the  entire  dielectric  loss  curve  in  the  frequency  domain  can  be  broken  up  into 
three  limiting  regions.  Near  the  loss  peak.  Debye  loss  mechanisms  are  still  taking  place  in  the  dipolar 
media,  and  so  the  classical  Debye  function  is  still  applicable.  However,  unlike  the  symmetric  toss  peak  of 
the  Debye  model,  asymmetric  power  law  tails  now  appear  in  the  loss  spectrum  both  at  higher  and  lower 
frequencies.  The  slopes  of  these  power  law  tails  depend  on  the  parameters  m  and  n.  In  this  model,  the 
relaxation  time  constant,  x,  is  defined  as  the  inverse  of  the  frequency  of  maximum  loss  and  is  determined 
from  the  parameters,  m.  n,  and  ay  The  quantity  coponly  equals  l it  when  n~l.  The  quantity  t  will  be  used 
in  temperature  scaling  of  the  structural  relaxation  in  the  data  analysis  below.  By  Fourier  transforming 
Equation  (11).  one  can  obtain  an  expression  for  the  depolarization  current  in  the  time  domain  in  terms  of 


the  confluent  hypergeometric  function 


42 


/(/)  =  *  P'(0)  exp(- copt)r" ,  Fx (l  - rn :  2  - » ; 


r(l  +  m-n) 
)r(l  +  m ) 


/(/)  cc  t  .  t«Mcop 

/(f) oc  exp(  co p f ),  t  =  \/a>p 
/(/)ocf'""\  t»\!  (Op 
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where  /  is  the  time  rate  of  change  of  P.  the  macroscopic  polarization  of  the  media.  Again,  there  is  a  Debye 
region  at  times  near  l/wp,  and  two  power  law  tails  at  times  shorter  and  longer  than  l/oy  The  temperature 

dependence  is  implicit  in  the  temperature  dependence  of  op. 

In  Fitiure  2.  the  KWW  and  D-H  models  are  compared  by  choosing  an  appropriate  log-log  plot.  It  is 
apparent  that  both  can  fit  the  same  data  depending  on  the  time  scale  examined.  In  addition,  if  the  power 
law  regions  have  nearly  equal  exponents  (symmetric  case),  then  the  two  models  may  be  difficult  to 
distinguish.  We  find  with  our  frequency  domain  data  below,  that  an  asymmetric  function  is  needed  to  fit 
the  data.  Thus,  the  D-H  model  is  more  general  than  the  KWW  model  and  is  able  to  fit  all  of  the  data. 

The  advantage  of  the  Dissado-Hill  model  over  the  previously  discussed  models  is  that  it  possesses  a 
stronger  physical  basis  being  consistent  with  causality  and  by  appropriately  accounting  for  temperature  and 
frequency  dependence  of  the  amplitude  of  the  response.  It  is  universal  in  that  it  has  been  applied  to  a  wide 
variety  of  dipolar  materials.39  Also,  it  provides  a  theoretical  description  for  dipolar  relaxation  in  both  the 
time  and  frequency  domains.  We  have  adapted  it  here  to  describe  both  dielectric  and  chielectnc  isothermal 

response  spectra. 


2.2.  Structural  Relaxation 

Models  for  structural  relaxation  account  for  the  temperature  dependence  of  the  characteristic  time  x 
found  in  the  reorientation  models  above.  Though  other  parameters,  namely  n  and  m.  are  necessary  to 
describe  the  time  dependence  of  the  reorientation,  they  are  not  as  strongly  temperature  dependent  as  t,  and 
no  good  models  exist  for  their  temperature  dependence.  As  above,  the  models  used  are  drawn  from  the 
glass  relaxation  literature  used  to  describe,  for  example,  dielectric  and  mechanical  relaxation.  There  have 
been  several  theoretical  models  introduced  to  describe  the  temperature  dependence  of  the  relaxation  time 
constant  based  on  entropy,  free  volume,  or  other  phenomenological  considerations.43  Of  these  models,  the 
three  most  widely  used  (along  with  closely  related  variations)  are.  the  Williams-Landel-Ferry  (WLF) 
model,  an  Arrhenius  model,  and  the  Adams-Gibbs  model.  The  last  model  is  the  most  general  in  that  it 
yields  the  other  two  in  appropriate  temperature  regions. 
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Stahelin  a/.44  investigating  polvetherimides  and  BPDA-PDA  polyimides  doped  with  lophtne 
chromophores.  have  found  that  x  is  correlated  to  the  glass  transition  temperature.  Ts  of  the  host  polymer  tn 
the  following  manner. 


/ 

r  =  r0  exp  ■ 
v 


B 

T  -T 

£o  L 


\ 

) 


(13) 


where  B  and  are  fitting  parameters,  and  T,  .T.+SO-C.  It  should  be  noted  that  this  type  of  temperature 
dependence  i,  similar  to  Wilhams-Landel-Ferry  (WLF)  o,  Vogel-Fulche,  (VF)  expressions  used  to 
describe  viscoelastic  properties  of  polymers  above  T,  ,  and  is  appropriate  only  in  that  temperature 

45. 4G 

range. 

Well  below  Tg,  the  nonlinear  optical  relaxation  was  found  to  follow  Arrhenius  behavior,  consistent 
with  other  relaxation  processes,  in  other  words. 


where  A  is  the  activation  energy  of  the  processes,  and  x0  is  a  parameter  describing  the  activation  entropy. 

In  the  Adam-Gibbs  entropic  model  of  the  glass  transition,  the  structural  relaxation  of  a  polymer  ts  the 
result  of  an  increase  in  the  number  of  molecular  segments  cooperatively  rearranging  as  the  temperature 

,  17.48.49 

decreases. 

Tool  originally  modeled  structural  relaxation  of  polymers  in  terms  of  Tf,  the  Active  temperature.  The 
method  for  calculating  Tf  for  a  given  polymeric  glass  was  introduced  by  Narayanaswamy5'  and  further 
developed  by  researchers  such  as  Hodge52  and  Moynihan.53'54  In  terms  of  Tf,  Hodge  described  the  Adam- 
Gibbs  type  temperature  dependence  of  the  relaxation  time  constant  as 


f 

T  =  A  exp 

V 


(15) 


or  scalinu  x  to  it’s  value  at  the  glass  transition  temperature. 
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where  t,  is  the  relaxation  time  constant  at  the  glass  transition  temperature.  B  is  an  activation  energy,  A  is  a 
time  scaling  parameter,  and  T2  is  the  zero  point  configurational  entropy  temperature.  The  physical  meaning 
of  the  fictive  temperature  is  shown  in  Figure  3.  At  temperatures  above  the  glass  transition  temperature,  Tf 

~T,  and  the  expression  for  x  reduces  to. 


=exp 


B 

V 


T-T, 


(17) 


which  is  just  the  Fulcher-Tammann-Hesse  or  the  equivalent  Williams-Landel-Ferry  expression.  Below  T, 
Tf  approaches  Tg  and 


r/rx  =  exp 


B 


Tg-T 


T 


(18) 


which  is  an  equivalent  expression  for  Arrhenius  temperature  behavior.  Thus,  the  Adam-Gibbs  model 
provides  a  more  general  description  of  the  temperature  dependence  oft.  because  it  can  be  used  to  describe 
relaxation  both  above  and  below  the  glass  transition  temperature. 


3.  Experimental  technique 

3. 1.  Material  Preparation 

The  polymer  systems  studied  in  this  work  included  two  guest-host  systems,  namely  DR1/PMMA  and 
DCM/PMMA,  and  one  side  chain  system.  DR1-MMA.56  The  glass  transition  temperatures  of  these 
materials  were  measured  using  differential  scanning  calorimetry  (DSC)  and  are  listed  in  Figure  4  along 
with  their  molecular  structures.  Samples  were  prepared  by  spin  coating  solutions  of  the  chromophore  and 
polymer  onto  ITO  coated  sapphire  substrates  and  baking  at  Tg+  30  °C  for  several  hours  yielding  2pm  thick 
films  that  consisted  of  10%  by  weight  of  chromophore.  Next,  two  coated  substrates  were  placed  in  a 
sandwich  configuration  and  fused  under  pressure  for  two  hours  at  a  temperature  above  the  polymer’s  glass 
transition  temperature.  The  pressure  rig  was  designed  to  apply  the  same  pressure  to  all  samples.  Then  the 
pressure  was  removed  and  the  sandwiched  film  was  annealed  for  one  hour  at  an  elevated  temperature  of 
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Tg+30C  to  relieve  any  built  up  stresses  in  the  sample.  Consistent  thermal  histories  were  applied  in  all 
cases. 

3.2.  Chromophore  reorientation  measurements 

Measuring  chromophore  orientational  relaxation  is  accomplished  by  monitoring  the  macroscopic 
polarization  of  a  NLO  polymer  system,  or  to  be  more  precise,  by  monitoring  the  linear. '/  } ,  and  nonlinear, 
X(2) ,  susceptibilities.  There  are  basically  four  types  of  experiments  that  can  be  performed  on  polymeric 
media  that  measure  these  tensors,  two  of  which  probe  x(l) ,  and  two  that  probe  X('’-  The  two  experiments 
that  probe  x(»  are  dielectric  loss  measurements  in  the  frequency  domain,  and  depolarization  current 
measurements  in  the  time  domain,  which  monitor  the  time  rate  of  change  ot  the  polarization.  As  for  the 
Xi2)  experiments,  these  are  SH  decay  measurements  in  the  time  domain  and  chielectric  relaxation 
measurements  in  the  frequency  domain.57  Time  domain  measurements  are  more  appropriate  at  long  time 
scales,  and  frequency  at  shorter  time  scales,  and  the  combination  allows  a  wide  range  of  time  scales  to  be 
probed.  To  minimize  sample  to  sample  measurement  variation,  all  four  experiments  were  performed  on  the 
same  sample.  To  analyze  the  experimental  results,  the  Dissado-Hill  model  was  chosen,  because  it  provided 
the  best  lit  to  both  the  time  domain  and  frequency  domain  data  in  all  cases. 

The  experimental  setup  used  for  time  domain  measurements  is  shown  in  Figure  5.  In  this 
configuration,  both  the  SH  light  generated  from  the  chromophores  and  the  depolarization  current  can  be 
measured  simultaneously.  SH  measurements  were  made  by  focusing  intense  laser  light  from  a  TuSapphtre 
source  onto  a  contact  poled  sample.  The  100  MHz  mode-locked  TiiSapphire  Argon  ion  pumped  laser 
produced  50  fs  pulses  as  measured  by  an  autocorrelator  with  an  average  output  power  of  approximately 
700mw  and  a  780nm  peak  wavelength.  After  the  fundamental  light  was  filtered  out,  the  SH  light  generated 
from  the  sample  was  collected  and  focused  onto  a  Hammamatsu  avalanche  photodiode.  The  electrical 
signal  from  the  photodiode  was  monitored  using  a  Stanford  Research  SRS850  dual  phase  lock  in  amplifier 
whose  reference  frequency  was  generated  from  an  optical  chopper.  The  depolarization  current  was 
monitored  using  a  Hewlett-Packard  6517  electrometer  that  also  served  as  the  voltage  source  for  poling  the 
sample. 


13 


Dielectric  measurements  were  made  using  a  CGA-83  raticarm  transformer  bridge.  The  frequency 
range  mtestigated  was  from  10  Hz  to  100  kHz  and  the  temperature  range  was  again  near  Tg  for  all 

samples. 

Chielectric  relaxation,  or  frequency  domain  electric  field  induced  second  harmonic  generation,  is  a 
relative l>  new  technique  introduced  by  J.A.  Cline  and  W.N.  Herman  to  observe  chromophore  relaxation  in 
the  frequency  domain.57  Their  method  is  similar  to  dielectric  relaxation  measurements  that  use  a  sinusoidal 
electric  Held  to  probe  the  linear  dielectric  response,  but  differs  in  that  it  measures  the  second  harmonic 
(SH)  light  generated  from  the  relaxing  chromophores  to  probe  the  nonlinear  optical  susceptibility.  The 
advantage  of  the  chielectric  relaxation  measurements  is  that  the  measured  SH  response  is  entirely  due  to  the 
chromophore.  while  in  the  dielectric  case,  both  the  chromophore  and  polymer  contribute  to  the  dielectric 
response.  It  is  also  easily  connected  to  the  time  domain  measurements  of  SHG  so  that  time  scales  over 
many  decades  can  be  probed. 

Using  a  mode-locked  Ti:Sapphire  laser,  a  more  quantitative  analysis  of  frequency  domain 
chromophore  relaxation  has  been  made.  Broadening  the  frequency  range  to  include  frequencies  from  10 
mHz  to  40  kHz.  has  allowed  measurements  to  be  made  at  temperatures  from  slightly  below  Tg  to 

temperatures  far  above  Tg. 

The  chielectric  frequency  domain  experimental  setup  is  nearly  the  same  as  the  time  domain  setup, 
except  for  a  few  variations  as  shown  in  Figure  6.  SH  measurements  were  made  by  applying  a  sinusoidal 
voltage  (  -400  V,,,,)  to  the  sample's  1TO  electrodes  while  the  intense  laser  light  from  a  TkSapphire  source 
was  focused  onto  the  sample.  After  the  fundamental  light  was  filtered  out.  the  SH  light  generated  from  the 
sample  uas  again  collected  and  focused  onto  a  Hammamatsu  avalanche  photodiode.  The  electrical  signal 
from  the  photodiode  was  monitored  using  a  Stanford  Research  SRS850  dual  phase  lock  in  amplifier  which 
also  served  as  the  sine  wave  generator  for  the  applied  voltage.  The  frequency  range  examined  in  this  work 
was  from  10  mHz  to  40  kHz.  and  the  temperature  range  was  near  Tg  for  all  samples. 

The  assumption  is  made  that  the  second  order  nonlinear  optical  susceptibility  varies  linearly  with 
respect  to  the  applied  electric  field,  which  means  if  the  applied  field  has  a  sinusoidal  time  dependence,  *(2) 
can  be  rewritten  as. 

=  A  sin(<y  /  +  <j>) 
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where  co  is  the  frequency  of  the  applied  field  and  (j)  is  the  phase  lag  of  the  chromophore  orientation  with 
respect  to  the  applied  field.  Since  the  SH  intensity  varies  as  the  square  of  x  > 

I2a  oc  (a  sin(<z>  t  +  ^))‘  =  ^(l  -  2cos(2(«  t  +  ^)))  •  (20) 

Thus,  by  setting  a  lock-in  amplifier  to  lock  in  at  twice  the  applied  voltage  frequency,  both  the  in  phase 

component  and  the  out  of  phase  component  of  can  be  measured  directly. 

The  fictive  temperature  was  determined  using  a  TA  Instruments  DSC  2910  differential  scanning 
calorimeter.  Rate  cooling  heat  capacity  measurements  were  performed  on  all  three  of  the  polymer  systems 
in  this  study  following  the  procedure  outlined  by  DeBolt  et  al.54  In  rate  cooling  heat  capacity  experiments, 
the  polymer  sample  is  first  heated  above  the  glass  transition  temperature,  then  cooled  below  Tg  at  a 
constant  rate,  and  finally  heated  at  20  °C/min  through  Tg.  This  thermal  cycle  was  repeated  for  several 
different  cooling  rates  ranging  from  20  °C/min  to  .5  °C/min.  From  the  heat  capacity  data  obtained  from 
these  sets  of  experiments,  the  temperature  dependence  of  the  fictive  temperature,  Tf,  was  determined  using 
the  relationship  between  the  heat  capacity  and  the  temperature  rate  of  change  of  Tf, 

dr,  c,(r)-c„(r) 

dT  Cp,(T)-CjT) 

where  CPE  and  Cpe  are  the  extrapolated  heat  capacities  in  the  glassy  and  equilibrium  temperature  regions. 

Results  and  Discussion 

Upon  removal  of  the  poling  field,  the  magnitude  of  x<2)  will  begin  to  decay.  This  decay  can  be 
monitored  by  measuring  the  intensity  of  SH  light  generated  from  the  relaxing  chromophores.  For  the  dye 
molecules  in  our  studies  there  are  only  two  nonzero  tensor  elements  of  y}2):  x(2)  3iiand  x(2>  333-  For  our 
measurements  we  arranged  the  geometric  setup  such  that  we  are  coupling  to  the  x(2)  333  element,  and  so  we 
have  the  following  expression  for  the  intensity  of  SH  light, 
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where  |-  '  is  the  SH  light  intensity.  p333  is  the  molecular  second  order  hyperpolarizability  tensor  element, 
and  9  is  the  angle  the  relaxing  dye  molecule  makes  with  the  direction  of  the  poling  field.38'59  The  final 
expression  for  XU)  hj  was  made  by  expanding  in  terms  of  Legendre  polynomials  and  then  approximating  in 
terms  of  <cos  0>  .‘,0  Note  that  Equation  (22)  indicates  that  the  dielectric  and  nonlinear  optic  experiments 
mainly  probe  the  same  order  parameter,  however,  the  nonlinear  optic  measurement  probes  only  the 
nonlinear  optical  chromophore.  An  example  decay  curve  for  DR1XPMMA  is  shown  in  Figure  7.  It  is  clear 
from  the  data,  that  neither  a  single  nor  double  exponential  would  fit  the  data.  At  the  same  time  the  SH 
decay  was  being  measured,  the  depolarization  current  caused  by  the  rotating  dipoles  of  the  relaxing 
chromophores  was  also  being  monitored.  In  theory  the  current  being  measured  should  be  entirely  due  to 
the  time  rate  of  change  of  the  macroscopic  polarization,  however,  low  frequency  dispersion  (LFD)  was 
found  in  all  of  the  samples  studied.  The  primary  source  of  LFD  is  probably  due  to  trapped  charges  being 
released.  An  example  of  the  depolarization  current  data  is  shown  in  Figure  8.  The  kink  at  long  times 
corresponds  to  the  LFD  due  to  trapped  charges. 

As  previously  mentioned,  dielectric  measurements  were  carried  out  in  the  frequency  range  from 
10  Hz  to  100  kHz.  A  detailed  account  of  the  basic  theory  behind  dielectric  relaxation  can  be  found  m  texts 
by  Jonscher39  or  McCrum.61  The  in  phase  and  out  of  phase  components  of  the  dielectric  spectra  for  the 
three  M.O  polymer  systems  studied  are  shown  in  Figure  9.  Both  a  and  P  peaks  are  visible  in  the  loss 
spectra.  The  peak  of  the  imaginary  part  of  the  dielectric  constant  determines  the  characteristic  frequency 
which  is  a  fit  parameter  in  the  D-H  model.  It  is  clear  that  the  characteristic  frequency  depends  strongly  on 
temperature.  This  dependence  will  be  compared  to  the  Adams-Gibbs  relaxation  model  below. 

Figure  10  shows  the  measured  X'3'  components  for  the  three  NLO  polymer  systems  in  this  study. 
We  see  close  agreement  between  the  dielectric  measurements  and  chielectric  spectra,  and  with  the  D-H 
model  fits.  We  also  determined  that  the  spectra  are  linear  in  the  electric  field  as  expected  from 
Equation  (19)  and  shown  in  Figure  11.  This  linearity  indicates  that  trapped  volume  or  surface  charges  are 
not  conmbuting  in  these  measurements.  In  time  domain  SHG  measurements,  the  trapped  charges  may 
affect  the  decay.  However,  we  have  carefully  shorted  the  electrodes  as  the  poling  field  was  turned  off  and 
did  not  see  any  obvious  effect  of  trapped  charges.  Further  evidence  for  this  is  given  below  where  the  time 
domain  and  frequency  domain  measurements  are  found  to  closely  correspond. 


16 


The  data  collected  from  all  four  sets  of  experiments  were  analyzed  within  the  framework  of  the 
Dissado- Hill  many  body  universal  relaxation  model.  As  previously  stated,  this  model  was  chosen  not  only 
because  it  provided  a  better  fit  to  experimental  data,  but  also  because  it  is  not  ent.rely  empirical,  and  has 
analytical  relaxation  expressions  for  both  the  time  and  frequency  domains.  In  calculating  the  relaxation 
time  constant,  t,  for  the  time  domain  measurements,  the  experimental  data  was  fit  to  the  time  domain 
Dissado-Hill  depolarization  current  express, on  from  Equat.on  (12).  From  this  the  three  Dissado-Hill 
parameters  n,  m,  and  cop  were  determined  and  inserted  into  Equation  (11)  to  calculate  the  frequency  of 
maximum  loss.  i.e.  the  reciprocal  of  t.  For  the  frequency  domain  experiments,  the  data  was  fit  directly  to 
Equation  (1 1)  to  calculate  x.  Putting  all  four  sets  of  experimental  data  together,  a  plot  of  the  temperature 
dependence  of  x  is  shown  in  Figure  12  for  the  DR1NPMMA  of  polymer  system.  It  should  be  noted  that 
while  the  dielectric  measurements  show  both  an  a  and  P  relaxation,  the  chielectric  data  shows  only  one 
relaxation  that  overlaps  the  dielectric  a  relaxation  data.  This  suggests  that  chromophore  relaxation  is 
strong  I  \  coupled  to  the  a  relaxation  of  the  polymer  host.  To  monitor  the  structural  relaxation  of  the 
polymer  host  the  temperature  dependence  of  the  fictive  temperature  is  calculated  from  the  rate  cooling 
enthalpy  experiments.  Example  enthalpy  results  are  shown  in  Figure  la. 

In  Figure  14.  the  temperature  dependence  of  x  for  all  three  polymer  systems  are  superimposed  on 
one  plot.  For  each  system,  the  temperature  is  scaled  to  the  respective  glass  transition  temperature.  By 
scaling  the  temperature  with  respect  to  Tg  of  each  system,  the  measured  x  values  from  the  a  relaxations 
appear  to  lie  on  a  single  master  curve  as  suggested  by  other  studies. 53  This  leads  to  predictive  capabilities 
for  the  useful  lifetime  of  the  electro-optic  properties  of  polymer-chromophore  systems.  It  is  also  apparent 
in  Figure  14  that  all  four  techniques  yield  consistent  results,  and  are  thus  likely  probing  the  same  physical 
mechanisms.  This  suggests  that  real  charge  effects  are  not  seen  under  these  experimental  conditions. 

Also  in  Figure  14.  the  Adam-Gibbs  theoretical  fit  with  Tf  calculated  from  the  DR1-MMA  enthalpy 
data  ,  along  with  Arrhenius  and  WLF  fits  have  been  superimposed  on  the  experimental  data.  The  WLF 
curve  fits  the  data  above  Tg  and  the  Arrhenius  curve  fits  the  data  below  Tg,  the  Adam-Gibbs  curve  fits  the 
data  in  all  temperature  ranges. 
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4.  Conclusion 

By  combining  results  from  both  time  and  frequency  domain  experiments,  a  complete  description  of 
chromophore  relaxation  in  polymeric  media  can  be  determined.  Using  the  Dissado-Hill  many  body 
universal  relaxation  model  to  extract  the  relaxation  time  constant,  t,  from  the  data,  and  the  Adam-Gtbbs 
formalism  to  describe  the  temperature  dependence  of  x,  we  have  found  good  agreement  between 
experiment  and  theory  over  13  orders  of  magnitude  in  time  and  over  temperatures  that  ranged  from  far 
below  T,  to  far  above  it.  We  have  also  shown  that  chromophore  relaxation  is  coupled  primarily  to  the  a 
relaxation  of  the  host  polymer.  In  addition,  after  scaling  the  temperature  to  Ts  and  the  relaxation  time 
constant,  r.  to  ts,  the  calculated  values  for  x  for  all  three  types  of  polymeric  systems  in  this  study,  seem  to 
have  the  same  reduced  temperature  dependence  which  indicates  this  type  of  relaxation  behavior  may  be 
universal,  at  least  within  the  same  types  of  glass  forming  polymers.  We  have  also  further  developed  the 
chielectric  relaxation  experimental  technique  for  measuring  chromophore  reorientation  in  the  frequency 
domain,  and  have  shown  how  it  can  be  used  to  directly  study  the  relaxation  mechanism  in  the 

chromophores. 

Currently  we  are  investigating  a  wider  variety  of  NLO  polymeric  media,  including  higher  Tg 
guest/host  systems,  polyimides.  and  also  cross-linked  polymers,  to  determine  if  the  results  found  in  this 
limited  >tudy  can  be  applied  in  general.  If  they  can.  the  techniques  discussed  here  may  prove  to  be  a  useful 
tool  in  predicting  which  NLO  media  will  be  viable  in  device  applications. 
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List  of  Figures 


Figure  1.  potential  energy  diagram  of  the  Dissado-Hill  many-body  two-level 
system.  The  two  preferred  orientation  states  are  represented  by  the  potential  energy 
wells.  The  shaded  regions  represent  energy  states  of  width  It,  resulting  from  particle 
interactions,  while  the  splitting  of  the  two  well  bottoms  is  due  to  the  applied  electric  field. 
Large  thermal  transitions  without  and  with  tunneling  assistance  over  the  energy  barrier  of 
height  A  are  represented  by  arrows  a  and  a'  respectively.  Arrow  b  represents  the  flip 
configurational  tunneling  transitions  while  arrow  c  represents  the  flip-flop  tunneling 
transitions. 

Figure  2.  Comparison  of  the  KWW  and  D-H  models  for  the  decay  of  the  polarization, 
P,  as  a  function  of  time.  At  short  times  the  two  models  are  nearly  identical,  however  at 
long  times  the  two  models  diverge. 

Figure  3.  q^e  temperature  dependence  of  the  specific  heat  of  a  glass  forming  liquid, 
where  ai  and  ag  represent  the  slopes  of  the  liquid  and  glassy  equilibrium  lines  which 
intersect  at  temperature  T2.  Cooling  the  sample  at  a  finite  cooling  rate,q,  the  fictive 
temperature,Tf,  can  be  graphically  determined  by  drawing  a  line  parallel  to  the  glassy 
equilibrium  line  from  the  particular  rate  cooling  curve  over  to  the  liquid  equilibrium  line. 
The  temperature  at  which  they  intersect  is  defined  as  the  fictive  temperature. 

Figure  4.  Molecular  structures  of  the  dyes  and  polymers  used  in  this  study,  and  their 
glass  transition  temperatures  as  calculated  by  differential  scanning  calorimetry. 

Figure  5.  Experimental  setup  for  time  domain  measurements. 


Figure  6.  Experimental  setup  for  frequency  domain  measurements. 

Figure  7.  Log-log  plot  of  isothermal  second  harmonic  decay  in  DR1\PMMA  at  50  °C. 
Figure  8.  Log-log  plot  of  isothermal  transient  current  decay  exhibiting  low  frequency 
dispersion  in  DR1\PMMA  at  50  °C. 

Figure  9.  Dielectric  in  phase  and  out  of  phase  spectra.  In  column  (a),  the  m  Phase 

dielectric  components  are  plotted  as  a  function  of  frequency  for  DCM\PMMA, 
DR1\PMMA,  and  DR1-MMA  respectively.  In  column  (b),  the  out  of  Phase  or  loss 
components  are  plotted  for  DCM\PMMA,  DR1YPMMA,  and  DR1-MMA  respectively.  In 
the  loss  spectra,  the  a  peaks  are  visibile  at  lower  frequencies  while  the  P  peaks  are  visible 
at  higher  frequencies. 

Figure  10.  chielectric  in  phase  and  out  of  phase  spectra.  In  column  (a),  the  in  Phase 

chielectric  components  are  plotted  as  a  function  of  frequency  for  DCM\PMMA, 
DR1\PMMA,  and  DR  1 -MM A  respectively.  In  column  (b),  the  out  of  phase  or  loss 
components  are  plotted  for  DCMVPMMA,  DR1\PMMA,  and  DR  1 -MM A  respectively. 
Figure  11.  Electric  field  dependence  of  chielectric  spectra  for  DCM\PMMA  at  fixed 
temperature.  In  (a)  and  (b)  the  in  phase  and  out  of  phase  components  of  y(2)  are  plotted 
as  a  function  of  frequency  for  various  poling  field  amplitudes.  In  (c)  and  (d)  the  in  phase 
and  out  of  phase  components  are  scaled  to  the  applied  poling  field  and  show  that  the  % 
response  is  linear  in  the  applied  field. 

Figure  12.  Reduced  temperature  dependence  plot  for  t  using  experimental  data  from 
transient  current,  second  harmonic,  dielectric,  and  chielectric  measurements  on 


DR1\PMMA. 


Figure  13.  Temperature  dependence  of  the  heat  capacity  for  DR1-MMA  after  rate 
cooling  at  5  °C/min.  The  thin  line  represents  the  extrapolated  heat  capacity  of  the  glassy 
region.  Cpg  while  the  thick  line  represents  the  extrapolated  heat  capacity  of  the 
equilibrium  region,  Cpe. 

Figure  14.  Reduced  temperature  dependence  plot  for  x  scaled  to  tg  for  the  three 
polymer  systems.  The  solid  line  represents  the  Adam-Gibbs  model  based  on  enthalpy 
measurements  from  DR1-MMA.  The  dashed  line  represents  an  Arrhenius  fit  based  on 
data  points  below  the  glass  transition,  while  the  dotted  line  represents  Fulcher-Tammann- 
Hesse  or  equivalently  WLF  behavior  above  the  glass  transition. 
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ABSTRACT 


We  have  investigated  crosslinkable  polyimides  for  both  passive  and  ac¬ 
tive  electro-optic  devices.  These  fluorinated  polyimides  are  soluble  in  the 
imidized  form  and  are  both  thermally  and  photo-crosslinkable  leading  to 
easy  processability  into  waveguide  structures  and  the  possibility  of  stable 
electro-optic  properties.  We  have  fabricated  channel  and  slab  waveguides 
and  investigated  the  mechanism  of  optical  propagation  loss  using  photother- 
mal  deflection  spectroscopy  and  waveguide  loss  spectroscopy,  and  found  the 
losses  to  arise  from  residual  absorption  due  to  the  formation  of  charge  trans¬ 
fer  states.  The  absorption  is  inhibited  by  fluorination  yielding  propagation 
losses  as  low  as  0.4  dB/cm  in  the  near  infrared.  Channel  waveguides  formed 
by  a  simple  wet  etch  process  are  observed  to  have  no  excess  loss  over  slab 
structures.  We  have  produced  electro-optic  polymers  by  doping  with  the 
nonlinear  optical  chromophores,  DCM  and  DADC;  and  a  process  of  con¬ 
current  poling  and  thermal  crosslinking.  Multilayer  structures  have  been 
investigated  and  poling  fields  optimized  in  the  active  layer  by  doping  the 
cladding  with  an  anti-static  agent. 


1.  INTRODUCTION 


Organic  polymers  have  received  much  attention  as  promising  materials 
for  integrated  optics  because  they  have  low  optical  losses  and  are  compati¬ 
ble  with  existing  photolithographic  fabrication  procedures.  [1)  They  can  be 


formed  into  both  passive  and  active  device  structures  and  are  easily  inte¬ 
grated  with  electronic  circuits  since  they  can  be  easily  spin  deposited  onto 
semiconductor  substrates.  Furthermore,  polymers  can  be  cast  into  multi¬ 
layers  containing  channels,  ribs,  and  other  routing  devices.  In  addition, 
for  active  applications,  their  low  dielectric  constant  implies  applicability  to 
microwave  frequency  electro-optic  devices.  [2,3]  Much  effort  in  improving 
electro-optic  polymers  is  currently  being  focused  on  optimizing  the  molec¬ 
ular  optical  nonlinearity,  the  fi/3  product,  and  in  increasing  the  thermal 
stability  of  the  electro-optic  coefficient. 

The  instability  of  the  electro-optic  coefficient  arises  from  the  electric 
field  poling  process  used  to  induce  the  polar  order  necessary  for  attaining  the 
Pockels  effect.  The  induced  polar  order  is  not  in  thermodynamic  equilibrium 
due  to  the  glassy  nature  of  the  polymer,  and  thus  orientational  relaxation 
occurs.  Current  approaches  to  optimizing  stability  has  included  the  use  of 
main  chain  chromophores,  side  chain  chromophores,  and  crosslinked  sys¬ 
tems.  [4-6]  These  materials  have  higher  glass  transition  temperatures  and 
reduced  orientational  mobility  as  compared  to  their  guest-host  analogs  and 
consequently  reduced  chromophore  relaxation.  Guest-host  polyimides  have 
been  employed  to  optimize  thermal  stability  of  induced  orientation  while 
maintaining  functionality.  Polyimides  are  well  known  in  the  electronics  in¬ 
dustry  for  their  thermal  integrity,  resistance  to  organic  solvents,  and  com¬ 
patibility  with  silicon  fabrication  processes.  Initial  guest-host  polyimides 
required  imidization  and  densification  at  high  temperature  requiring  that 
dopant  molecules  have  high  thermal  decomposition  temperatures.  [7,8]  In 
studies  where  the  chromophores  did  not  thermally  decompose,  sublimation 
and  plasticization  occurred.  [9,10]  Alternatives  to  high  temperature  imidiza¬ 
tion  have  been  developed  and  include  chemical  imidization  and  preimidized 
soluble  polyimides.  [11,12]  In  these  cases,  densification  for  extended  peri¬ 
ods  of  time  was  necessary  to  minimize  the  free  volume  around  the  chro¬ 
mophore  and  reduce  thermal  relaxation  of  induced  orientation.  [13]  More 
recently,  side-chain  polyimides  have  been  developed.  [14-16]  These  materi¬ 
als  imidize  at  temperatures  below  the  thermal  decomposition  temperatures 
of  most  chromophores  and  have  demonstrated  excellent  thermal  stability  of 
induced  orientation  for  over  200  hours  at  elevated  temperatures.  The  ad¬ 
ditional  orientational  stability  is  presumably  due  to  the  hindered  rotational 
motion  of  the  large  polyimide  chains. 

We  review  here  our  approach  to  the  development  of  device-quality 
electro-optic  materials  using  perfluorinated  preimidized  fully  aromatic  poly¬ 
imides.  Preiraidization  makes  processing  doped  polymer  systems  more  flex¬ 
ible,  by  making  functionalizing  easier  while  allowing  spin  coating  of  soluble 
fully-imidized  polymers.  Fluorination  increases  solubility  while  decreasing 
optical  loss  and  refractive  index.  Fully  aromatic  polymers  allow  for  the 


best  high  temperature  properties.  We  found  that  by  introducing  alkylated 
aromatic  crosslinking  groups  which  can  be  photo-  or  thermally  activated, 
optical  losses  were  greatly  reduced,  and  high  quality  waveguides  could  be 
fabricated.  The  photo-crosslinking  also  permits  a  simple  liquid  etch  process 
to  define  waveguides,  and  provides  a  chemical  hook  to  which  chroraophores 
may  be  covalently  attached  to  these  polymers.  We  present  our  efforts  to 
incorporate  this  material  system  in  microwave  frequency  integrated  optic 
devices. 


2.  POLYMER  WAVEGUIDES 


Waveguiding  in  passive  polymer  structures  were  studied.  Single-layer 
guides  on  glass  or  Si02  substrates  were  studied  to  determine  the  intrinsic 
optical  properties  such  as  refractive  index  and  loss.  Polymer  core-cladding 
structures  were  also  fabricated  in  order  to  evaluate  materials  in  passive  device 
geometries. 

We  studied  optical  losses  in  polyimide  waveguides  in  order  to  under¬ 
stand  the  mechanisms  responsible  for  losses  and  to  develop  new  polyimides 
with  lower  losses.  [17,18]  Low  loss  waveguides  could  be  produced  by  attach¬ 
ing  perfluoro  (CF3)  groups  to  the  polymer  backbone.  The  addition  of  these 
groups  decreases  absorption  in  the  visible  wavelengths  by  reducing  charge 
transfer  and  creating  polymer  deformation. 

Polyimides  were  synthesized  with  increasing  amounts  of  fluorination 
as  shown  in  Fig.  1  (in  order  of  increasing  fluorination).  The  details  of  syn¬ 
thesis  have  been  reported  elsewhere.  [19]  Refractive  indices  and  optical  losses 
were  measured  by  waveguide  techniques  as  well  as  by  photothermal  deflec¬ 
tion  spectroscopy.  [17,18]  Refractive  indices  were  observed  to  decrease  with 
increasing  fluorination.  The  low  refractive  index  results  from  its  low  den¬ 
sity  due  to  the  presence  of  alkylated  aromatic  crosslinking  groups  [20]  and 
a  high  degree  of  fluorination.  It  was  determined  that  propagation  losses 
in  waveguides  was  mainly  due  to  absorption  arising  from  charge  transfer 
absorption  bands  in  the  polyimide.  Losses  were  observed  to  decrease  with 
increasing  fluorination.  Because  of  the  low  losses  and  fabrication  advantages 
of  cross-linkable  materials,  Ultradel  9020D  and  9120D  were  selected  for  fur¬ 
ther  studies  of  both  active  and  passive  waveguides.  Table  1  contains  physical 
and  optical  properties  of  Ultradel  9020D  and  9120D.  It  is  seen  from  the  ta¬ 
ble  that  these  materials  are  indeed  promising  materials  for  both  passive  and 
active  applications.  They  exhibit  both  low  dielectric  and  optical  losses,  are 
thermally  stable,  and  possess  useful  coefficients  of  thermal  expansion.  In  ad- 
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FIG.  1.  Chemical  structures  of  polyimides  used  in  this  study  are  listed 
in  order  of  increasing  fluorination.  Ultradel  9020D  and  9120D  are  crosslink- 
able  polyimides  that  can  be  used  as  a  core  cladding  pair.  (Ri  =  alkylated 
aromatic  crosslinking  group  as  described  in  reference  18) 
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Table  1.  Physical  and  optical  properties  of  Ultradel  9020D  and  9120D. 
(TGA=thermal  gravimetric  analysis:  CTE=coefficient  of  thermal  expansion) 


dition,  the  refractive  indices  suggest  that  a  core-clad  structures  with  9120D 
core  and  9020D  cladding  would  be  possible  to  make. 

Clad  passive  waveguide  structures  were  fabricated  using  a  wet  etch 
process  which  we  will  describe  in  a  later  section.  Structures,  as  shown  Fig.  2. 
with  Ultradel  9120D  as  a  core  and  9020D  as  a  cladding  were  made  and 
waveguide  losses  were  measured.  Losses  were  also  observed  to  be  highly 
cure  dependent  as  shown  in  that  figure.  At  temperatures  greater  than  350° 
C.  increased  losses  are  due  to  polyimide  decomposition.  However,  as  seen 
in  the  figure,  losses  in  the  core-clad  waveguides  as  low  a  0.4  dB/cm  were 
attained.  In  addition,  these  losses  were  as  low  as  observed  in  single  layer 
waveguides,  indicating  that  the  waveguide  fabrication  process  introduced  no 
additional  losses.  [17] 

Nonlinear  optically  active  materials  are  made  by  mixing  nonlinear  op¬ 
tical  chromophores  into  the  host  material.  The  introduction  of  nonlinear 
chroraophores  may  introduce  additional  losses  in  the  form  of  scattering  sites 
and  absorption  tails  which  can  be  further  enhanced  upon  poling.  [21]  Fig.  3 
depicts  the  structure  of  our  nonlinear  optical  dopants,  DCM  [22,23]  and 
DADC  [24].  In  both  cases  significant  losses  are  introduced  by  the  chro- 
mophore.  Fig.  4  depicts  the  additional  waveguide  losses  associated  with 
incorporation  of  the  nonlinear  optical  dye.  This  loss  arises  from  absorption 
tails  induced  in  the  materials,  probably  due  to  the  formation  of  charge  trans¬ 
fer  complexes  between  the  dye  and  polymer.  The  fully  aromatic  nature  of 
the  polymer  likely  contributes  to  this.  The  formation  of  these  absorption 
tails  may  be  inhibited  by  using  structures  less  likely  to  form  charge  transfer 
complexes. 


3.  ELECTRO-OPTIC  PROPERTIES 


Electro-optic  polymers  are  produced  by  doping  the  polyimide  with 
nonlinear  linear  optical  chromophores  and  followed  by  electric  field  poling. 
[25]  Ultradel  9020D  was  doped  with  thermally  stable  chromophores,  DCM 
and  DADC,  and  electric  field-poled  to  impart  an  electro-optic  response. 
DCM  is  a  highly  nonlinear  chromophore  with  low  molecular  weight  that 
is  photo- bleach  able  in  the  UV.  [22]  Using  DCM  as  a  dopant,  we  can  ob¬ 
tain  large  nonlinearities  because  high  number  densities  are  possible.  While 
DCM’s  decomposition  temperature  remains  above  the  processing  tempera¬ 
ture  of  Ultradel  9020D,  other  researchers  have  reported  sublimation  at  tem¬ 
peratures  as  low  as  220°  C.  [10]  DADC,  a  chromophore  developed  for  high 
temperature  electro-optic  applications,  has  a  high  nonlinearity,  high  molec- 
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FIG.  2.  Core-clad  waveguide  structure.  Waveguide  loss  data  at  A  =  800 
nm  plotted  as  a  function  of  curing  temperature.  Thermal  crosslinking  occurs 
near  200°C. 
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FIG.  3.  Chemical  structures  and  properties  of  thermally  stable  high  tem¬ 
perature  chromophores  DADC  and  DCM. 
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FIG.  4.  PDS  spectra  showing  the  absorptive  losses  for  Ultradel  9020D, 
Ultradel  9120D,  DCM/Ultradel  9020D,  and  DADC/Ultradel  9020D.  Chro- 
mophore  loading  is  17%  by  weight  fraction.  Curing  temperature  for  these 
samples  were  300  C,  300C,  175  C  and  300  C,  respectively 
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Table  2.  Dependence  of  the  electro-optic  coefficient  on  poling  temperature. 
The  maximum  electro-optic  coefficient  for  DCM  and  DADC  immediately 
after  poling  is  scaled  to  unity. 


ular  weight,  and  is  also  photo-bleachable  in  the  UV.  [24]  Using  DADC  we 
expect  lower  nonlinearities  than  DCM  because  of  lower  number  densities 
at  a  given  weight  fraction.  D ADC’s  larger  structure  eliminates  sublimation 
even  at  temperatures  as  high  as  400°  C.  TGA  (thermogravimetric  analysis) 
experiments  have  shown  only  a  2%  weight  loss  in  nitrogen  at  temperatures 
as  high  as  400°  C.  [22] 

The  optimum  poling  temperature  was  determined  by  comparing  sam¬ 
ples  poled  with  the  same  field  at  varying  temperatures.  Orientation  and 
relaxation  were  probed  by  measuring  the  in-plane  electro-optic  coefficient 
at  800  nra  with  a  Mach-Zehnder  interferometer.  [26]  The  electro-optic  co¬ 
efficient  was  measured  immediately  after  poling  for  both  sets  of  samples. 
Afterwards,  the  samples  were  placed  into  an  oven  at  125°  C  for  approxi¬ 
mately  50  hours,  and  the  electro-optic  coefficient  was  remeasured.  Table  2 
shows  the  results  for  both  DADC  and  DCM.  For  DADC  the  electro-optic 
coefficient  before  aging  is  constant  with  temperature.  After  thermal  aging, 
the  sample  poled  below  200°  C  decayed  while  those  above  200°  C  remained 
stable.  This  indicates  that  crosslinking  occurs  just  above  200°  C. 

Table  2  shows  that  the  ratio.  r33(t)/r33(0),  for  DCM  is  fairly  constant 
for  temperatures  above  200°  C.  At  higher  temperatures,  ^(O)  was  lower.  We 
attribute  this  to  sublimation  of  the  chromophore  at  elevated  temperatures. 
Interestingly,  the  chromophore  retention  depends  on  the  initial  concentration 
which  may  be  evidence  of  the  chromophore  plasticizing  the  polymer.  The 
enhanced  orientational  decay  of  DCM  at  all  temperatures  may  be  attributed 
to  plasticizing  of  the  polyimide.  [12] 

The  equation  relating  the  electro-optic  coefficient  to  its  microscopic 


parameters  at  a  particular  wavelength  is  given  by  [25], 
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where  N  is  the  chromophore  number  density,  Pzzz  is  the  second  order  mi¬ 
croscopic  susceptibility,  L^(p)  is  the  third-order  Langevin  function,  fi  is  the 
dipole  moment,  Ep  is  the  poling  field,  kg  is  the  Boltzmann  constant,  Tp  is 
the  poling  temperature.  f°  is  the  Onsager  local  field  factor  for  the  static 
poling  field.  f  UJf2uJ,fUJ'  are  the  Lorenz-Lorentz  local  field  factors  at  optical 
frequencies,  uf  and  2 uJ  are  the  frequencies  used  in  the  second  harmonic 
generation  experiment,  uj0  and  uj  are  the  resonance  and  electro-optic  exper¬ 
iment  frequencies,  respectively.  Fig.  5  shows  the  agreement  between  Eq.  1 
and  experiment  indicating  full  poling  and  verifies  the  electronic  nature  of  the 
electro-optic  effect.  The  solid  lines  define  the  (+/— )15%  error  in  the  mea¬ 
surement  of  fi(3.  The  slight  offset  of  electro-optic  coefficient  at  larger  poling 
fields  is  probably  due  to  increased  birefringence  with  increasing  poling  fields. 


4.  THREE-LAYER  ACTIVE  WAVEGUIDES 


Nonlinear  optical  materials  must  be  formed  into  channel,  slab,  or  rib 
waveguides  for  integrated  optical  applications.  The  resistance  of  crosslinked 
Ultradel  9020D  to  organic  solvents  makes  it  a  suitable  material  for  layered 
devices.  We  formed  three-layer  slab  waveguides  by  spinning  and  curing  suc¬ 
cessive  layers  of  Ultradel  9020D  on  ITO-coated  glass.  The  layered  structures 
consisted  of  DCM-  and  DADC-doped  Ultradel  9020D  cores  surrounded  by 
neat  Ultradel  9020D  cladding  layers.  To  obtain  noncentrosymmetric  or¬ 
dering  as  well  as  large  poling  field  induced  nonlinearities  in  three-layered 
devices,  it  is  essential  for  the  cladding  layers  to  be  more  conductive  than  the 
core  layer  so  that  the  voltage  drop,  and  hence  the  poling  field,  across  the 
active  layer  is  as  large  as  possible.  Direct  measurements  of  the  poling  field 
are  not  possible  in  three-layer  samples  because  the  voltage  splitting  between 
core  and  cladding  layers  is  unknown.  The  effective  poling  field  can  be  deter¬ 
mined  indirectly  by  measuring  the  macroscopic  susceptibility’s  dependence 
on  the  externally  applied  field  and  comparing  to  theory.  The  macroscopic 


FIG.  5.  Electric  field  dependence  of  electro-optic  coefficient  for  17% 
DADC  by  weight  fraction.  Solid  lines  represent  (+/-)15%  uncertainty  in 
fx/3  values  as  input  to  the  theory.  Data  points  are  from  electro-optic  experi¬ 
ment  where  r33  =  3ri3  was  assumed. 

susceptibility  was  measured  using  the  Rotational  Maker  Fringe  (RMF)  tech¬ 
nique  as  a  function  of  poling  field  for  single-  and  three-layer  samples.  [27] 
Gold  electrodes  were  sputtered  onto  the  samples  for  contact  poling.  After 
contact  poling,  the  gold  was  removed  and  the  RMF  experiment  was  per¬ 
formed.  Measurements  were  made  using  the  2nd  Stokes  line  (A  =  1.367/x)  of 
a  hydrogen-filled  Raman  cell  that  was  pumped  by  the  output  of  a  Nd:YAG 
pumped  pulsed  dye  laser.  The  measured  second  harmonic  coefficients  in  the 
single  layer  films  for  both  dopants  was  found  to  agree  with  an  oriented  gas 
model  similar  to  Eq.  1.  The  agreement  between  theory  and  experiment 
shows  that  the  oriented  gas  model  adequately  describes  the  poling  induced 
ordering  in  polyimide  materials. 

The  data  for  three-layer  samples  appear  in  Fig.  6  for  DCM-  and 
DADC-doped  samples.  The  two  sets  of  data  points  on  the  graphs  repre¬ 
sent  different  methods  of  calculating  the  poling  field  across  the  active  layer. 
The  open  data  points  assume  the  resistivities  of  all  three  layers  are  identical 
and  the  resulting  poling  field  is  obtained  by  dividing  the  applied  voltage  by 
the  total  thickness.  The  shaded  data  points  use  a  scaling  factor  to  fit  the 
experimental  data  to  the  theoretical  prediction  of  Eq.  3  with  the  poling 
field  as  an  adjustable  parameter.  If  we  assume  the  three-layer  sample  can 
be  modeled  by  resistors  in  series,  then  the  scaling  factor  can  be  related  to 
the  cladding  and  core  resistances  by  the  following  equation, 


V*f,  = 


Rr 


Rcore  F  2  Rcladding 


-vn 


applied 


(3) 


where  RCOre  and  Rdaddmg  are  the  resistances  of  the  core  and  cladding  layers, 
and  Veff  and  Vapplied  are  the  effective  and  applied  poling  voltages.  [28]  The 
open  data  points  show  the  inadequacy  of  assuming  equal  resistivities  to 
determine  the  poling  field  across  the  active  layer.  For  DCM  triple-layers,  a 
value  of  U=0.4  is  used  for  the  best  fit  adjustable  parameter,  while  for  DADC 
triple-layers,  a  value  of  0.11  is  obtained.  Conductivity  of  thin  single-layer 
films  was  used  to  independently  calculate  the  effective  poling  field.  The 
conductivity  and  thickness  of  DADC  triple-layers  poled  at  250°  C  predicts 
a  correction  factor  of  0.1  in  excellent  agreement  with  the  best  fit.  Thus, 
the  voltage  division  model  appears  valid.  For  DCM  triple-layers  poled  at 
225°  C  Eq.  3  predicts  an  effective  poling  field  an  order  of  magnitude  smaller 
than  the  best  fit  adjustable  parameter.  This  discrepancy  may  be  due  to 
migration  of  DCM  into  the  cladding  layers  which  lowers  the  conductivity  of 
the  core  layer  while  increasing  the  conductivity  of  the  cladding  layers.  This 
migration  is  consistent  with  the  sublimation  observed  in  single-layer  films. 
The  optimization  of  poling  requires  that  the  conductivity  of  the  cladding  be 
much  larger  than  the  core  to  maximize  the  electric  field  in  the  active  layer. 

To  this  end,  we  have  incorporated  ionic  dopant  [29]  consisting  of  a  graft 
polymer  incorporating  a  quaternary  ammonium  salt  (used  as  an  anti-static 
agent)  into  the  cladding  layers  to  increase  the  conductivity.  A  polymeric  salt 
distributes  throughout  the  cladding  layers  in  contrast  to  small  molecule  salt 
which  would  reside  at  the  surface.  Measurements  of  the  imaginary  part  of 
the  dielectric  constant  which  is  proportional  to  the  conductivity  are  shown 
in  Fig.  7  as  a  function  of  concentration.  It  is  seen  at  low  concentrations 
that  the  conductivity  is  nearly  linear,  but  becomes  exponential  at  higher 
concentrations.  When  the  axis  is  scaled  to  the  measured  conductivities  of 
the  undoped  polvimide  and  the  DADC  doped  polyimides.  the  conductivity 
of  the  core  layer  is  shown  as  a  dotted  line.  We  expect  that  the  ionic  doped 
layer  at  18%  loading  will  have  the  requisite  conductivity  to  amplify  the  field 
in  the  core  substantially.  Fig.  8  depicts  the  second  harmonic  coefficient  for  a 
three  layer  film  consisting  of  18%  ionic  doped  cladding  and  a  core  containing 
15%  DADC  by  weight.  The  line  is  a  fit  to  the  oriented  gas  model.  The  fit 
factor  is  about  0.3  indicating  that  the  field  is  divided  nearly  equally  between 
the  three  layers.  This  factor  is  considerably  less  than  the  value  of  0. 5-0.6 
we  would  expect  from  the  conductivity  in  Fig.  7  though  substantially  higher 
than  the  0.11  without  the  ionic  dopant.  This  discrepancy  may  be  due  to  the 
decomposition  or  sublimation  of  the  ionic  dopant  during  the  poling  process. 
Refractive  index  measurements  of  the  doped  polymer  before  and  after  high 
temperature  processing  support  this  conclusion. 


FIG.  6.  <^33  versus  poling  field  for  three  layer  stacks  of  (a)  DCM/Ultradel 
9020D  and  (b)  DADC/Ultradel  9020D.  Solid  lines  represents  fit  to  an  ori¬ 
ented  gas  model.  Data  points  are  measurements  from  three- layer  samples 
using  the  poling  field  as  an  adjustable  parameter.  Best  results  (shaded  data 
points)  are  obtained  using  0.4  as  the  best  fit  parameter  for  DCM  samples 
and  0.1  for  DADC  samples,  (unshaded  data  points)  assume  that  resistances 
of  all  layers  are  equal. 
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FIG.  7.  Imaginary  dielectric  constant  for  ionic  doped  polyimide.  The 
dotted  line  indicates  the  approximate  core  equivalent  conductivity. 
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FIG.  8.  Three  layer  second  harmonic  coefficient  data  for  three  layer  films 
with  18%  by  weight  ionic  doping  in  the  core  and  15%  DADC  doping  in  the 
core.  Fit  is  to  an  oriented  gas  model. 
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FIG.  9.  Processing  schematic  for  fabrication  of  single  mode  polyiraide 
waveguides. 


5.  CHANNEL  DEVICE  FABRICATION 


A  processing  scheme  for  creating  channel  devices  has  been  developed. 
This  process  is  shown  in  Fig.  9.  The  photosensitive  polymer  is  spin-deposited 
onto  a  substrate,  and  soft-cured  (below  200°C)  to  remove  solvent.  The  poly¬ 
mer  is  exposed  through  a  photomask  leading  to  photo-crosslinking,  and  then 
partially  etched  with  organic  solvent  developers  into  a  rib  or  channel  pattern. 
The  coating  is  then  post-baked  at  300  —  350°C  to  remove  residual  solvent 
and  to  fully  crosslink.  Either  rib  or  channel  waveguides  can  be  produced 
with  this  method. 

Fig.  10  shows  a  SEM  micrograph  of  a  rib  structure  fabricated  using 
photosensitive  polyimides  showing  their  high  quality  structure.  Channels 
and  ribs  4-5  microns  wide  have  been  produced.  After  post-cure,  the  films 
can  withstand  processes  such  as  metallization,  overcoating,  thermal  cycling, 
and  solvent  exposure.  Waveguide  loss  spectroscopy  was  carried  out  on  three- 
layer  rib  waveguides  composed  of  two  variations  of  Ultradel  9020D  polymers 
with  appropriate  refractive  indices.  No  excess  optical  loss  was  introduced  by 
channel  waveguide  fabrication. 


FIG.  10.  SEM  of  a  5.0  /zra  y-branch  waveguide  etched  in  Ultradel  9020D. 


6.  RF  ELECTRODE  DESIGN 


We  describe  here  electrode  design  for  our  polyimide  materials  similar 
to  that  carried  out  for  other  polymer  systems.  [2,30]  In  the  design  of  high 
speed  modulators,  a  key  factor  is  the  relative  velocities  of  the  microwave  and 
optical  signals  as  they  propagate  down  the  device.  The  coupling  efficiency 
between  the  waves  depends  on  the  velocity  match  between  the  waves  with  the 
upper  frequency  limit  depending  on  this  phase  velocity  match.  Polyimides 
are  well-known  for  their  low  dielectric  constant;  in  fact,  the  microwave  di¬ 
electric  constant  and  the  optical  refractive  index  squared  can  be  very  close 
so  that  the  velocity  match  is  not  the  limiting  factor  in  such  devices.  We 
have  designed  raicrostrip  electrodes  for  propagation  of  microwave  frequency 
signals  in  a  modulator.  [31] 

The  geometry  of  the  microstrip  line  and  the  optical  waveguide  are 
shown  in  Figs.  11  and  12  with  the  microstrip  line  placed  directly  over  the 
polymer  waveguide.  An  impedance  match  taper  is  also  included.  To  char¬ 
acterize  the  structure,  the  microstrip  characteristic  impedance  Zo,  the  mi¬ 
crostrip  propagation  constant,  7  =  0  +  j/?,  the  length  of  the  coupling  region 
L,  the  source  voltage  Vo,  and  the  source  and  load  impedances  Z5  and  Zi 
are  required.  In  the  case  of  this  geometry,  the  total  optical  phase  shift  due 
to  the  electro-optic  effect  is  given  by, 


M  -  6r(f)  sin (0  -  fmt) 


(4) 


(5) 


where  6  is  the  electro-optic  phase  shift  given  by 

~n2rV0TL 

s  =  ^jr~ 

with  n  being  the  index  of  refraction,  r  the  electro-optic  coefficient,  H  the 
device  thickness,  T  the  electrical-optical  overlap  parameter,  A  the  optical 
wavelength,  and  L  the  propagation  length.  The  parameter  <p  is  the  phase 
delay  between  the  RF  and  optical  wave  which  will  be  defined  below.  The 
frequency  dependent  amplitude  response  is  given  by, 

R(f)  =  10  log  r(f)  -1-  Prf  (6) 

where  Prf  is  the  required  RF  power,  and 
/  r \  __  /  a2  +  b2 

”  V  e2aL(l  -  p5pLe~2ctL)2{{K L)2  +  (aL)2)2 

Here,  we  have  ps  and  pi  the  RF  source  and  load  reflection  coefficients, 
K  =  (27r/ An)  /c,  with  An  =  y/^eff  —  n  proportional  to  the  difference  in 
phase  velocities  between  the  microwave  and  optical  wave.  Finally, 

a  =  L(K  cos(KL)  +  asin (KL)  -  KeaL) 

+  LT(e2aLK  cos  (KL)  -  e2ocLasm(KL)  —  K)  (8) 

b  =  L(K  sin  (KL)  —  acos(KL)  +  aeaL) 

+  LT(e2aLK  sm(KL )  +  e2aL  a  cos(K  L)  —  a)  (9) 

and 

T  =  pLpse~2aL.  (10) 


The  phase  delay  is  (b  =  tan-1  (a/6).  Because  An  is  so  small  in  poly- 
imides,  it  the  attenuations  and  reflections  that  limit  the  3  dB  bandwidth. 
The  responsivity  R(f)  was  calculated  for  a  number  of  heights,  H,  as  a  func¬ 
tion  of  frequency  and  shown  in  Fig.  13.  In  generating  Fig.  13,  a  50  ohm 
impedance  at  the  input  and  output  ports,  a  width  of  15  /xm,  a  length  of  1.2 
cm,  dielectric  constant  =  2.6,  tan  6  —  0.002,  and  refractive  index  of  1.61  were 
used.  Note  the  3  dB  bandwidth  out  to  50  GHz.  The  microstrip  fabricated 
using  these  parameters  predicts  an  R(f)  electro-optic  response  as  shown  in 
Fig.  14.  These  predictions  were  made  from  measured  insertion  losses.  Good 
agreement  is  seen  with  the  simulations,  and  a  bandwidth  greater  than  40 
GHz  is  expected. 
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FIG.  11.  Section  of  microstrip-coplanar  transition  on  polyiraide  waveguide. 


FIG.  12.  Cross  section  of  the  optical  waveguide  raicrostrip  device. 


Frequency  (GHz) 

FIG.  13.  Calculated  modulator  amplitude  response  versus  frequency  for 
a  microstrip  width  of  15  /xra. 

7.  CONCLUSIONS 


We  have  described  the  development  of  a  crosslinkable  polyimide  material 
system  for  integrated  optical  applications.  First,  we  optimized  the  opti¬ 
cal  and  waveguide  processing  properties,  and  later  introduced  functionality 
by  guest-host  inclusion  and  electric  field  poling.  From  optical  waveguide 
and  PDS  measurements  we  determined  the  primary  mechanism  for  optical 
waveguide  loss  in  polyimide  waveguides  is  absorption  from  electronic  absorp¬ 
tion  tails.  The  measured  waveguide  loss  in  Ultradel  9020D  of  0.4  dB/cm  at 
800  nm  and  predicted  losses  of  0.4  dB/cm  at  1300  nm  make  these  materials 
excellent  candidates  for  integrated  optical  applications.  Introduction  of  the 
active  chromophores,  DCM  and  DADC,  increased  the  optical  losses  as  deter¬ 
mined  by  PDS.  We  believe  that  these  long  absorption  tails  may  be  a  general 
problem  in  many  electro-optic  polymers,  but  especially  so  in  aromatic  ones 
like  those  studied  here. 

We  then  investigated  the  electric  field  poling  conditions.  The  optimum 
poling  temperature  was  found  to  be  slightly  below  the  sublimation  temper¬ 
ature  of  DCM  and  well  below  the  thermal  decomposition  temperature  of 
most  thermally  stable  chromophores.  Both  guest-host  systems  display  ex¬ 
cellent  thermal  stability  of  the  electric  field-induced  orientation  at  room 
temperature.  At  elevated  temperatures,  significant  orientational  relaxation 
occurred  in  DCM-doped  samples.  The  enhanced  decay  may  be  attributed 
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FIG.  14.  Responsivity,  R(f)  of  a  polymer  electro-optic  modulator. 


to  plasticizing  of  the  polyimide  and  greater  rotational  mobility  of  DCM. 
The  DADC-doped  polyimide  retained  75%  of  its  initial  orientation  after  200 
hours  at  125°  C. 

Three-layer  samples  of  doped  polyimide  showed  nonlinearities  smaller 
than  comparable  single-layer  films.  The  decrease  in  nonlinearity  is  attributed 
to  lower  poling  fields  due  to  voltage  division  across  the  three-layer  sample  and 
possibly  chromophore  migration  into  cladding  regions.  The  voltage  splitting 
can  be  eliminated  by  substituting  existing  cladding  layers  with  more  con¬ 
ductive  cladding  layers  as  we  observed  with  an  ionic  dopant  in  the  cladding 
layers. 

We  have  also  described  a  simple  waveguide  processing  procedure  for 
forming  channel  and  rib  structures  which  takes  advantage  of  the  pho¬ 
tocrosslinking  nature  of  polyimide  materials  and  does  not  introduce  excess 
optical  losses.  We  have  also  described  the  fabrication  of  modulating  elec¬ 
trodes  with  a  predicted  frequency  response  to  at  least  40GHz.  We  believe 
that  this  material  system  is  an  excellent  candidate  for  active  and  passive 
integrated  optical  applications.  Farther  refinement  in  processing  and  refrac¬ 
tive  index  control  are  all  that  is  necessary  for  passive  applications.  For  active 
applications,  the  identification  of  suitable  chromophores  which  do  not  intro¬ 
duce  absorption  is  necessary.  Covalent  attachment  of  active  chromophores 
to  crosslinking  sites  should  be  straightforward,  and  would  lead  to  vastly 
improved  performance  regarding  optical  nonlinearity  and  stability. 
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